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Zoonotic  microbes  have  historically  been,  and  continue  to  emerge  as,  threats  to  human  health.  The  recent 
outbreaks  of  highly  pathogenic  avian  influenza  virus  in  bird  populations  and  the  appearance  of  some  human 
infections  have  increased  the  concern  of  a  possible  new  influenza  pandemic,  which  highlights  the  need  for 
broad-spectrum  detection  methods  for  rapidly  identifying  the  spread  or  outbreak  of  all  variants  of  avian 
influenza  virus.  In  this  study,  we  demonstrate  that  high-density  resequencing  pathogen  microarrays  (RPM) 
can  be  such  a  tool.  The  results  from  37  influenza  virus  isolates  show  that  the  RPM  platform  is  an  effective 
means  for  detecting  and  subtyping  influenza  virus,  while  simultaneously  providing  sequence  information  for 
strain  resolution,  pathogenicity,  and  drug  resistance  without  additional  analysis.  This  study  establishes  that 
the  RPM  platform  is  a  broad-spectrum  pathogen  detection  and  surveillance  tool  for  monitoring  the  circulation 
of  prevalent  influenza  viruses  in  the  poultry  industry  and  in  wild  birds  or  incidental  exposures  and  infections 
in  humans. 


Recent  outbreaks  of  Nipah  virus,  severe  acute  respiratory 
syndrome  virus,  and  avian  influenza  virus  reiterate  the  impor¬ 
tance  of  zoonotic  microbes  as  potential  threats  to  human 
health  (26).  Influenza  virus  causes  particular  concern,  owing  to 
the  repeated  nature  of  influenza  pandemics  and  their  potential 
to  result  in  significant  mortality,  exemplified  by  the  1918  influ¬ 
enza  pandemic.  To  date,  most  influenza  A  virus  subtypes  (e.g., 
H2N2  and  H10N7)  resulting  from  combinations  of  the  16  hem¬ 
agglutinin  (HA)  and  9  neuraminidase  (NA)  types  are  detected 
in  wild  birds  and  poultry  (9,  40).  In  a  survey  of  38,609  samples 
from  wild  birds,  55  different  HA  and  NA  subtype  combinations 
were  identified,  with  H4N6  appearing  as  the  most  prevalent 
subtype,  followed  by  H7N7  and  H6N2  (24).  The  emergence  of 
H5N1  since  1997  in  Asia,  the  Middle  East,  Europe,  and  Africa 
amplifies  concerns  about  the  broad  natural  diversity  of  host 
species  (mostly  aquatic  and  migratory  birds)  which  provide 
rapid  geographical  distribution  of  new  strains  and  enable 
transmission  to  human  populations  (24,  25).  Recent  major 
outbreaks  in  domestic  poultry  and  wildfowl  populations  caused 
by  different  serotypes,  including  H5N1,  H5N2,  H7N1,  H7N3, 
H7N4,  and  H7N7,  indicate  that  the  threat  is  not  from  a  single 
serotype  (25). 
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Detection  and  discrimination  of  all  potential  influenza  A 
virus  subtypes  is  needed  to  identify  the  introduction  of  zoo¬ 
notic  strains  to  humans,  monitor  the  status  of  these  pathogens 
in  their  natural  hosts,  and  minimize  epidemic  spread  if  trans¬ 
missible  human  infections  occur.  An  effective  surveillance 
assay  could  rapidly  detect  and  identify  all  subtypes  of  avian 
influenza  virus  and  provide  useful  secondary  information  re¬ 
lated  to  specific  functional  mutations  which  alter  pathogenicity 
or  drug  resistance.  For  example,  the  low-pathogenicity  H5N1 
viruses  should  be  differentiated  from  the  highly  pathogenic 
strains  by  a  mutation  in  an  HA  cleavage  site  (a  multibasic 
cleavage  motif,  PQRERRRKKRG),  a  deletion  of  20  amino 
acids  in  the  NA  protein,  and  a  signature  amino  acid  substitu¬ 
tion,  E627K,  in  the  PB2  protein  (1,  23). 

Viral  culture  paired  with  serological  HA  typing  is  the  current 
standard  method  for  detecting  and  typing  influenza  A  viruses. 
These  procedures  are  time-consuming,  taking  days  or  even 
weeks  to  provide  specific  results.  Several  molecular  diagnostic 
approaches  including  reverse  transcription  (RT)-PCR,  real¬ 
time  PCR,  PCR-enzyme-linked  immunosorbent  assay,  and 
spotted  oligonucleotide  microarrays  provide  fast  and  sensitive 
alternatives  to  viral  culture  (5,  8,  12-15,  20,  23,  27,  31,  33,  34, 
38).  While  promising,  these  methods  either  are  limited  to  de¬ 
tecting  only  a  few  subtypes  or  provide  a  very  limited  range  of 
genetic  resolution.  Additional  time-consuming  characteriza¬ 
tion,  such  as  direct  sequencing,  is  required  for  analysis  of  strain 
variations  and  specific  mutations  that  contribute  to  or  predict 
influenza  virus  pathogenicity,  host  range,  drug  resistance,  and 
vaccine  efficacy. 

Alternative  strategies,  such  as  the  use  of  RT-PCR  coupled 
with  pyrosequencing  (7,  28),  RT-PCR-electrospray  ionization 
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(ESI)-mass  spectrometry  (MS)  (30),  or  resequencing  pathogen 
microarrays  (RPM)  (6, 16,  17, 19,  36),  allow  tracking  of  genetic 
changes  and  supply  subspecies  identification.  The  pyrose- 
quencing  technique  is  currently  limited  to  short  fragments  and 
is  applied  to  detection  of  H5N1  or  selected  drug  resistance 
markers  (3,  4,  7,  28).  The  RT-PCR-ESI-MS  method,  devel¬ 
oped  for  detection  of  all  avian  serotypes,  has,  to  date,  only 
demonstrated  tracking  of  genetic  changes  in  human  influenza 
virus  samples.  The  RPM  technology  is  the  only  one  of  these 
technologies  currently  under  development  for  simultaneous 
detection  and  identification  of  influenza  A  virus  variants  to¬ 
gether  with  a  large  number  of  other  viral  and  bacterial  patho¬ 
gens  that  may  elicit  similar  flulike  illnesses.  Furthermore,  the 
RPM  technology  separates  and  partially  decouples  the  ampli¬ 
fication  of  limiting  templates  by  multiplex  RT-PCR  from  the 
selection  of  microarray  contents  and  detection  capability, 
which  alleviates  constraints  on  primer  selection  while  still  pro¬ 
viding  the  required  specificity. 

Herein,  we  investigate  the  performance  of  new  versions  of 
the  respiratory  pathogen  microarray  (TessArray  RPM-Flu  3.0 
and  3.1,  subsequently  designated  RPM-Flu  (see  Table  SI  in 
the  supplemental  material)  for  detection  and  differential  iden¬ 
tification  of  all  subtypes  of  the  influenza  A  virus  HA  and  NA 
genes  in  a  single-pass  assay.  Previous  studies  demonstrated  the 
ability  of  RPM  technology  to  detect  targeted  pathogens  with 
analytical  and  clinical  sensitivities  and  specificities  that  are 
similar  to  (or  improved  over)  those  for  existing  technologies, 
while  simultaneously  providing  sequence  information  for 
strain  resolution  (6, 16, 17, 19,  22,  36).  The  RPM-Flu  arrays  are 
designed  and  constructed  to  allow  thorough  coverage  of  86 
bacterial  and  viral  agents,  including  respiratory  pathogens  and 
zoonotic  organisms  considered  to  be  significant  risks  for  hu¬ 
man  health,  e.g.,  severe  acute  respiratory  syndrome  virus. 

About  30%  of  the  RPM-Flu  array  is  dedicated  to  targeting 
all  16  HA  and  9  NA  alleles  of  avian  influenza  A  viruses.  The 
HA  and  NA  genes  represented  on  the  microarray  are  based  on 
prevalent  strains  of  influenza  A  viruses  circulating  in  the  avian 
population  within  the  last  few  years.  Previous  studies  have 
shown  that  a  single  sequence  on  a  resequencing  microarray 
could  reliably  detect  and  serotype  strains  with  as  much  as  10  to 
15%  variation  (16,  18,  19,  36).  Thus,  this  technology  has  the 
potential  to  detect  all  possible  influenza  A  virus  HA  and  NA 
subtype  combinations.  The  data  presented  here  demonstrate 
that  the  RPM  platform  is  an  effective  means  for  universal 
detection  and  identification  of  all  subtypes  of  avian  influenza  A 
viruses  and  provides  useful  secondary  information  related  to 
pathogenicity  and  drug  resistance. 

(This  work  was  presented  in  part  at  the  107th  general  meet¬ 
ing  of  the  American  Society  for  Microbiology,  2007.) 

MATERIALS  AND  METHODS 

Influenza  A  virus  HA  and  NA  templates.  HA1,  HA3,  NA1,  and  NA2  sequences 
were  amplified  from  total  nucleic  acids  extracted  from  live,  trivalent,  nasally 
administered  influenza  vaccine  (FluMist  2004/05;  Medimmune,  Inc.,  Gaithers¬ 
burg,  MD).  The  influenza  A  virus  H5N1  (Asian  lineage)  real-time-PCR/RT- 
PCR-positive  control  (catalog  no.  VA2711;  Centers  for  Disease  Control  and 
Prevention  Laboratory  Response  Network)  was  kindly  provided  by  the  Centers 
for  Disease  Control  and  Prevention  (Atlanta,  GA).  The  influenza  A  virus  HA2, 
HA4-16,  and  NA3-9  templates  were  synthesized  by  BlueHeron  Biotechnology, 
Inc.  (Bothell,  WA)  (see  Table  S2  in  the  supplemental  material). 


AVIAN  INFLUENZA  VIRUS  IDENTIFICATION  989 

Specimen  collection  and  processing.  Cloacal  swabs  from  migratory  birds 
(mostly  waterfowl),  cloacal  and/or  tracheal  swabs  from  commercial  poultry,  tra¬ 
cheal  and  lung  tissue  samples  from  dead  birds,  human  throat  swabs,  and  one 
human  lung  sample  taken  from  a  deceased  patient  were  collected  by  Naval 
Medical  Research  Unit  3  (NAMRU3;  Cairo,  Egypt).  These  samples  were  cul¬ 
tured  using  chicken  eggs  and/or  MDCK  cells.  Sample  collection  and  viral  culture 
techniques  were  performed  as  described  in  the  World  Health  Organization 
Manual  On  Animal  Influenza  Diagnosis  and  Surveillance,  Version  2002.5,  revi¬ 
sion  1  (39).  Total  RNA  was  extracted  from  culture  isolates;  1:100  to  1:1,000 
dilutions  were  used  for  molecular  diagnosis. 

Microarray  design.  The  details  of  the  design  and  selection  process  for  the  gene 
targets  for  RPM-Flu  arrays  were  described  in  a  previous  publication  (37).  Briefly, 
after  target  genes  were  chosen  based  on  the  association  between  sequence 
diversity  and  clinically  relevant  phenotypic  diversity,  in  silico  modeling  (21)  was 
employed  to  select  a  set  of  probes  that  specifically  recognize  only  the  targeted 
clade  of  organisms  while  accounting  for  all  or  most  of  the  recognized  variants 
(i.e.,  strains  or  subtypes)  within  that  clade.  The  RPM-Flu  arrays  (TessArray 
RPM-Flu  3.0  and  3.1;  TessArae,  LLC,  Potomac  Falls,  VA)  were  designed  to 
maximize  coverage  of  respiratory  pathogens,  targeting  all  16  HA  and  9  NA  alleles 
of  avian  influenza  A  viruses  as  well  as  three  matrix  (M)  alleles,  one  nonstructural 
protein  ( NS1 )  allele,  and  one  polymerase  basic  subunit  2  ( PB2 )  region  (see  Table 
SI  in  the  supplemental  material).  About  30%  of  the  RPM-Flu  arrays  are  dedi¬ 
cated  to  the  detection  of  influenza  viruses.  Because  numerous  respiratory  patho¬ 
gens  needed  to  be  detected  using  a  limited  amount  of  microarray  space,  most  HA 
and  NA  genes  were  represented  by  a  single  sequence  optimized  for  prevalent 
strains  of  influenza  A  viruses  circulating  in  the  avian  population  within  the  last 
few  years  (see  Table  S2  in  the  supplemental  material).  We  have  focused  on 
ensuring  the  typing  of  the  most  common  variants,  while  ensuring  identification 
(with  potentially  less  detailed  typing  results)  for  less  common  variants.  There  are 
instances  where  less  common  types  will  have  more  than  15%  variation  for  some 
targets.  This  will  not  prevent  detection  of  the  influenza  virus  but  may  prevent 
subtyping  to  a  single  type  and  instead  result  in  a  list  of  several  equally  likely 
subtype  identifications.  Definite  subtyping  then  relies  on  other  targets  on  the 
array.  The  H10N7  samples  discussed  in  Results  serve  as  an  example  of  this.  The 
NA7  tile  is  designed  based  on  the  H7N7  subtype,  and  there  is  approximately  20% 
difference  in  the  NA  sequences  of  the  H7N7  and  the  H10N7  subtypes.  For  a 
sample  of  H10N7,  the  sequence  generated  from  the  NA7  tile  provides  detection 
but  indicates  both  H10N7  and  H7N7  as  being  equally  possible.  The  information 
from  the  HA10  (positive)  and  the  HA7  (negative)  subtypes  allow  a  final  deter¬ 
mination  that  H10N7  has  been  detected.  For  each  of  the  HA3,  HA6,  HA1,  NA1, 
and  NA8  genes,  two  different  reference  sequences  were  used  because  multiple 
sequence  alignments  revealed  two  major  prevalent  groups  that  could  not  be 
represented  by  a  single  reference  sequence. 

Primer  design.  Four  independent  multiplex  primer  mixes,  separated  to  sim¬ 
plify  primer  design  and  multiplex  optimization,  were  used  to  amplify  a  total  of 
187  targeted  sequences  represented  on  RPM-Flu  arrays  (298  primers).  The 
gene-specific  primer  pairs  for  all  targets  on  the  RPM-Flu  chips  (see  Table  SI  in 
the  supplemental  material)  were  designed  according  to  the  criteria  described 
previously  (16,  17).  To  accommodate  the  genetic  variation  of  RNA  viruses  and 
ensure  amplification  of  closely  related  variants,  a  software  script  automated  a 
primer  selection  algorithm  developed  by  our  group  (unpublished  data,  available 
upon  request)  to  select  primers  from  defined  primer  regions  of  each  target. 
Briefly,  target  sequences  were  used  to  search  against  GenBank  to  find  all  se¬ 
quences  sharing  at  least  80%  sequence  identity.  Consensus  sequences  that  cov¬ 
ered  the  primer  regions  flanking  the  target  sequences  were  generated  from  the 
sequences  found.  Primers  were  chosen  to  be  between  18  and  25  bp  in  length,  with 
predicted  annealing  temperatures  between  55  and  60°C  and  at  most  one  wobble 
nucleotide  (M,  R,  S,  W,  Y,  or  K).  When  necessary,  two  primer  pairs  instead  of 
one  were  used  to  amplify  one  target  sequence. 

Of  the  four  multiplex  primer  mixes,  one  was  dedicated  to  the  influenza  virus 
targets,  including  the  aforementioned  HA  and  NA  alleles  and  the  M  genes  from 
H3N2  and  H1N1,  which  served  as  general  markers  of  the  presence  of  influenza 
A  viruses.  The  M,  NS1,  and  PB2  genes  from  the  H5N1  subtype  were  included  to 
allow  more-detailed  analysis  of  critical  regions  subject  to  functional  mutation 
(see  Table  SI  in  the  supplemental  material). 

Multiplex  RT-PCR  amplification.  The  multiplex  RT-PCRs  were  performed  as 
previously  described  (16),  with  the  following  modifications.  For  the  RT  step, 
primer  LN  was  replaced  by  primer  NLN  (a  random  9-mer  with  a  linker  se¬ 
quence).  One  picogram  each  of  two  internal  controls  (NAC1  and  triosephos- 
phate  isomerase)  and  4  |xl  of  the  total  nucleic  acids  extracted  from  either  clinical 
specimens  or  laboratory  controls  were  used.  The  RT  reaction  products  were  split 
into  four  5-|xl  volumes  for  four  different  multiplex  PCRs.  Platinum  Taq  DNA 
polymerase  (Invitrogen  Life  Technologies,  Carlsbad,  CA)  was  replaced  by 
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GoTaq  DNA  polymerase  (Promega  Corporation,  Madison,  WI)  in  the  PCR. 
Linker  primer  NL  instead  of  primer  L  was  used  with  50  to  150  nM  each  of 
primers  from  mixes  1  to  4  in  the  multiplex  PCR.  The  amplification  reaction  was 
carried  out  with  an  initial  incubation  at  25°C  for  10  min  and  then  preliminary 
denaturation  at  94°C  for  2  min;  followed  by  16  cycles  at  94°C  for  30  s,  45  to  60°C 
for  30  s  (with  an  incremental  increase  of  1°C  per  cycle),  and  72°C  for  90  s;  and 
then  24  cycles  at  94°C  for  30  s  and  60°C  for  120  s.  The  amplified  products  from 
all  four  PCRs  were  combined  into  a  single  volume  and  subjected  to  purification 
and  processing  prior  to  hybridization  to  the  RPM-Flu  chips. 

Microarray  hybridization  and  analysis.  The  overall  procedure  of  microarray 
processing  protocol  was  based  on  a  previously  published  procedure  (see  Fig.  SI 
in  the  supplemental  material).  Microarray  hybridization  and  processing  and 
image  scanning  and  processing  for  production  of  FASTA  output  files  were 
performed  as  previously  described  (16).  Final  pathogen  identification  was  per¬ 
formed  using  Computer-Implemented  Biological  Sequence  Identifier  version  2.0 
software  (22),  an  automatic  pathogen  identification  algorithm  based  on  nucleic 
acid  sequence  alignment,  which  was  developed  and  tested  in  detail  in  previous 
studies  (16,  17). 

Reference  assays.  Hemagglutination  inhibition  assays  (HI)  were  performed  at 
NAMRU3,  following  the  World  Health  Organization  Manual  On  Animal  Influ¬ 
enza  Diagnosis  and  Surveillance,  version  2002.5,  revision  1  (39).  RT-PCR- 
ESI-MS  was  performed  as  previously  described  (30)  at  the  Naval  Health 
Research  Center  (NHRC;  San  Diego,  CA).  Additional  assays,  including 
conventional  and/or  real-time  RT-PCR,  were  also  performed  as  further  confir¬ 
mation  at  the  NHRC  and/or  the  Naval  Research  Laboratory  (NRL;  Washington, 
DC).  Real-time  RT-PCR  was  conducted  to  quantify  the  number  of  avian  influ¬ 
enza  virus  genomes  present  in  each  sample  (32).  Furthermore,  additional  sub- 
type-specific  RT-PCR  amplification  assays  were  performed  on  discordant  sam¬ 
ples  on  the  basis  of  both  RPM-Flu  and  RT-PCR-ESI-MS  results. 

Full-length  amplifications  of  M  and  NS1  genes,  using  the  universal  primer  set 
described  by  Hoffmann  et  al.  (11),  were  also  performed  on  discordant  samples  to 
obtain  unbiased  de  novo  sequence  results.  When  the  universal  primer  set  did  not 
produce  sufficient  material  for  conventional  sequencing,  subtype-specific  RT- 
PCR  amplifications  were  performed  based  on  both  RPM-Flu  and  RT-PCR- 
ESI-MS  results.  Amplified  products  from  HA,  M,  NA,  and  NS1  genes  were 
purified  and  sent  to  Macrogen  USA  (Gaithersburg,  MD)  for  Sanger/electro- 
phoresis-based  sequencing  using  either  correspondent  universal  or  specific 
primers. 

Nucleotide  sequence  accession  numbers.  All  nucleotide  sequences  reported 
in  this  study  are  available  at  GenBank  (http://www.ncbi.nlm.nih.gov/GenBank 
/index.html)  under  accession  numbers  EU599273  to  EU599315. 


RESULTS 

Analytical  sensitivity.  The  limited  representation  of  various 
HA  and  NA  subtypes  in  available  clinical  samples  constrained 
efforts  to  demonstrate  the  limits  of  detection  for  all  HA  and 
NA  alleles.  To  validate  the  RPM  platform  for  full-range  influ¬ 
enza  virus  detection,  artificially  generated  gene  fragments  (900 
to  1,730  bp)  (see  Table  S2  in  the  supplemental  material)  were 
used  as  an  alternate  means  to  estimate  platform  limits  of  de¬ 
tection.  The  analytical  sensitivities  for  unequivocal  detection 
and  differentiation  of  all  the  HA  1-16  and  AM  1-9  sequences 
represented  on  the  RPM-Flu  were  between  10,000  and  1,000 
genome  copies  per  assay  with  various  templates  (data  not 
shown).  This  sensitivity  range  with  test  templates  was  previ¬ 
ously  shown  to  correspond  to  excellent  detection  performance 
with  clinical  samples  (17).  These  results  indicated  that  the 
RPM-Flu  approach  should  be  an  effective  means  of  detecting 
and  subtyping  influenza  virus  directly  from  clinical  samples. 

Identification  of  influenza  A  virus  isolates.  Flaving  demon¬ 
strated  the  capability  of  RPM-Flu  for  detecting  and  differen¬ 
tiating  all  HA  and  NA  subtypes,  a  study  was  undertaken  to 
characterize  37  cultured  avian  influenza  virus  specimens  col¬ 
lected  in  the  Middle  East.  This  study  effectively  compared  the 
utility  of  the  RPM  platform  with  respect  to  the  traditional  HI 
for  avian  influenza  virus  detection.  The  samples  were  addition¬ 


ally  tested  using  alternate  PCR  methods  available  at  the  time 
(RT-PCR-ESI-MS  and  PCR  assays,  carried  out  at  the 
NHRC).  These  samples  were  blinded  with  respect  to  all  prior 
results  and  tested  at  the  NRL  by  using  the  RPM-Flu  platform. 
Sample  identities  were  revealed  only  after  the  RPM  results 
were  finalized  (Table  1). 

For  the  21  samples  identified  as  the  H5  serotype  by  HI,  all 
methods  produced  identical  results.  Three  additional  samples, 
for  which  the  HI  results  were  not  available,  were  H5N1  posi¬ 
tive  by  PCR.  Two  samples  identified  as  influenza  A  virus  by  HI 
were  negative  by  all  other  methods,  which  suggested  that  sam¬ 
ple  degradation  and/or  low  titers  may  have  been  involved.  Of 
the  remaining  11  samples,  6  samples  were  identified  as  other 
influenza  virus  serotypes  and  5  were  identified  only  as  influ¬ 
enza  A  virus  by  HI.  The  PCR  tests  (specifically  targeting  in¬ 
fluenza  A  virus,  H1N1,  H3N2,  and  H5N1)  provided  only  spe¬ 
cies  identification  for  these  samples.  Of  the  six  samples  for 
which  HI  made  serotype  identifications,  RT-PCR-ESI-MS 
identified  two  samples  as  H5N1,  three  samples  as  H7N7,  and 
one  sample  as  influenza  A  virus,  which  contradicted  the  H7, 
H9,  H10,  and  Hll  identifications  made  by  HI.  RPM-Flu  made 
the  same  serotype  identifications  as  HI  for  five  samples,  and 
one  sample  was  identified  only  as  influenza  A  virus  (Table  1). 
Of  the  five  samples  where  HI  offered  only  species  identifica¬ 
tion,  RT-PCR-ESI-MS  (four  samples)  and  RPM-Flu  (three 
samples)  offered  serotype  identification  or  identified  them  as 
influenza  A  virus.  For  the  two  samples  where  both  methods 
offered  serotype  identifications,  the  results  did  not  agree  (Ta¬ 
ble  1). 

Confirmation  of  influenza  virus  isolates  by  sequence  anal¬ 
ysis.  Conventional  sequencing  of  HA,  M,  NA,  and  NS  I  by  use 
of  universal  primers  (11)  was  attempted  on  the  11  discordant 
samples  and  2  H5Nl-positive  samples  as  controls.  When  the 
universal  amplification  procedure  failed,  subtype -specific  prim¬ 
ers  were  employed  based  on  the  typing  identifications  made  by 
both  the  RPM-Flu  and  the  RT-PCR-ESI-MS  platforms.  In  the 
event  that  conventional  sequencing  was  successful,  sequences 
generated  for  M  and  NS1  genes  did  not  provide  conclusive 
information  for  subtype  identification,  but  direct  HA  and  NA 
gene  sequencing  agreed  with  the  results  from  the  RPM-Flu 
and  HI  platforms.  In  cases  where  HI  failed  to  provide  conclu¬ 
sive  identification  but  the  RPM-Flu  provided  detection,  the 
identification  was  concordant  with  the  de  novo  sequencing 
analysis  (Table  1). 

Identification  of  highly  pathogenic  influenza  A  strains.  Trans¬ 
lation  of  RPM-Flu  resequencing  data  for  the  HA  and  PB2 
sequences  of  the  20  avian  H5Nl-positive  samples  confirmed 
that  13  samples  contained  the  HA  cleavage  site  and  15  samples 
indicated  key  amino  acid  substitutions  in  the  PB2  protein. 
Eight  samples  had  base  calls  before  and  after  the  location  of 
the  NA  deletion.  In  these  samples,  this  is  a  strong  indication 
that  a  deletion  is  present.  The  remaining  samples,  unfortu¬ 
nately,  had  base  calls  only  after  the  location  of  the  NA  deletion, 
so,  while  suggestive,  it  is  possible  that  the  absence  of  base  calls 
before  this  location  is  associated  with  the  low  titer  of  the 
sample  and  not  with  the  presence  of  a  deletion.  The  overall 
results  indicated  that  75%  of  the  H5N1  samples  had  at  least 
one  highly  pathogenic  strain  marker  and  no  conflicting  indica¬ 
tions.  The  remaining  samples  failed  to  provide  base  calls  in  all 
the  marker  locations;  therefore,  no  clear  determination  of  low- 
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TABLE  1.  Identification  of  influenza  A  virus  by  RPM,  HI,  and  RT-PCR-ESI-MS 


Sample  no. 

Sample  origin 

Result6  for: 

HI  (PI) 

RPM 

RT-PCR-ESI-MS 

PCR 

De  novo  sequencing 

2006900207 

Human 

H5 

H5N1 

H5N1 

H5N1 

— 

2006900588 

Poultry" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900755 

Human 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900756 

Human 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900839 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900840 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900843 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900844 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900845 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006901162 

Great  egret 

H5 

H5N1 

H5N1 

H5N1 

- 

2006902782 

Human 

H5 

H5N1 

H5N1 

H5N1 

- 

2006902786 

Human 

H5 

H5N1 

H5N1 

H5N1 

- 

2006902834 

Human 

H5 

H5N1 

H5N1 

H5N1 

- 

2006902838 

Human 

H5 

H5N1 

H5N1 

H5N1 

- 

2006902991 

Human 

- 

H5N1 

H5N1 

H5N? 

- 

2006902992 

Human 

- 

H5N1 

H5N1 

H5N? 

- 

2006903458 

Human 

- 

H5N1 

H5N1 

H5N1 

- 

2006906089 

Poultry" 

H5 

H5N1 

H5N1 

H5N? 

- 

2006906375 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006906608 

Duck" 

H5 

H5N1 

H5N1 

H5N1 

- 

2005909464 

Goose" 

H5 

H5N1 

H5N1 

H5N1 

- 

2005909467 

Duck" 

H5 

H5N1 

H5N1 

H5N1 

- 

2006900590 

Poultry" 

H5 

H5N1 

H5N1 

H5N1 

H?N1 

2006902764 

Chicken" 

H5 

H5N1 

H5N1 

H5N1 

H5N1 

2006905588 

Goose 

H7 

H7N7 

Flu  A 

H?N1 

H7N7 

2004909864 

Shoveler 

Flu  A 

H7N7 

H9N2 

Flu  A 

H7N? 

2004900600 

Shoveler 

H10 

H10N7 

H7N7 

Flu  A 

H10N7 

2005912823 

Teal 

Flu  A 

H10N7 

H7N7 

Flu  A 

H10N7 

2005912908 

Teal 

H10 

H10N7 

H7N7 

Flu  A 

H10N7 

2004900845 

Shoveler 

H10 

H10N7  or  H10N5 

H5N1 

Flu  A 

H10N7 

2004900688 

Teal 

Hll 

H11N? 

H7N7 

Flu  A 

H11N? 

2005912306 

Yellow-legged  bird 

Flu  A 

H13N6  or  H13N8 

NEW 

Flu  A 

H13N? 

2005909888 

Teal 

Flu  A 

Flu  A 

H5N2 

Flu  A 

- 

2003920431 

Teal 

H9 

Flu  A 

H5N1 

H?N2 

- 

2005910801 

Pintail 

Flu  A 

Flu  A 

H5N3 

Flu  A 

- 

2005910876 

Teal 

Flu  A 

NEG 

Low  amp  (H3N2) 

NEG 

- 

2006902947 

Human 

Flu  A 

NEG 

NEG 

Nl/NEG" 

- 

"  Domestic. 

b  NEG,  negative;  Flu  A,  influenza  A  virus;  Low  amp,  low-level  amplification;  — ,  no  PCR  results;  PI,  first  passage  after  culture. 
c  PCR  results  provided  by  the  NRL. 


or  high-pathogenicity  strains  could  be  made.  The  HA  peptide 
translated  from  the  two  H7N7  samples  indicated  low  virulence 
strains  (data  not  shown). 

Drug  resistance  information.  In  addition  to  tracking  the 
mutations  that  were  indicative  of  HPAI,  RPM  results  provided 
the  ability  to  track  drug  resistance  markers.  Previous  work  has 
shown  that  the  NA  protein  may  be  characterized  by  four 
critical  amino  acid  substitutions,  El  19V,  R152K,  H274Y,  and 
R292K,  associated  with  resistance  against  neuraminidase  in¬ 
hibitors.  The  amino  acid  substitution  R292K  was  of  particular 
significance,  as  this  mutation  results  in  complete  resistance  to 
oseltamivir  and  zanamavir  (2).  The  results  obtained  from 
translation  of  the  NA  gene  from  RPM-Flu  resequencing  anal¬ 
ysis  of  the  H5Nl-positive  samples  confirmed  that  21  avian 
H5N1  samples  with  good  base  call  rates  harbored  the  E119V 
mutation.  Furthermore,  the  R152K  mutation  was  not  observed 
in  17  out  of  21  samples,  the  R292K  mutation  in  20  out  of  21 
samples,  and  the  H274Y  mutation  in  9  out  of  21  samples  (data 
not  shown).  These  data  indicated  that  most  of  the  avian  H5N1 
samples  identified  were  presumably  sensitive  to  neuraminidase 


inhibitors.  Specific  single-amino-acid  substitutions  in  any  one 
of  four  critical  amino  acids  of  the  M2  protein,  L26F,  V27A/T, 
A30T/V,  and  S31N/R,  have  been  associated  with  amantadine 
resistance  (10,  29).  Unfortunately,  the  results  obtained  from 
translation  of  the  M2  genes  did  not  give  conclusive  results  for 
amantadine  resistance  (data  not  shown). 

DISCUSSION 

The  data  presented  herein  demonstrate  the  capabilities  of  a 
resequencing  microarray  approach  for  influenza  virus  surveil¬ 
lance.  The  combination  of  bioinformatics  (target  design  and 
pathogen  identification  algorithms)  and  amplification  strate¬ 
gies  applied  in  the  RPM  platform  maximizes  the  potential  for 
successful  detection  and  facilitates  molecular  epidemiological 
surveillance  by  identifying  and  discriminating  prevalent  sub- 
types.  The  accuracy  of  the  resequencing  microarrays  dem¬ 
onstrates  utility  not  only  for  detection  and  identification  of 
emerging  and  novel  subtypes  of  avian  influenza  viruses  but  also 
for  tracking  of  pathogenicity  and  drug  resistance  mutations.  In 
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light  of  current  concerns  over  the  potential  emergence  and 
spread  of  highly  pathogenic  avian  influenza  viruses  among 
humans  leading  to  a  pandemic,  the  RPM-Flu  assay  offers  an 
ideal  monitoring  method,  as  it  can  detect  all  of  the  prevalent 
serotypes  and  resolve  critical  mutations.  This  assay  potentially 
allows  for  specific  and  early  intervention  during  outbreaks, 
thereby  reducing  the  potential  economic  impact  associated 
with  persistent  disease.  Noteworthy  is  the  fact  that  reassort¬ 
ment  of  various  HA  and  NA  alleles,  which  can  readily  expand 
avian  influenza  virus  diversity  and  promote  the  emergence 
of  epidemiologically  important  new  strains,  can  be  detected 
readily  by  the  RPM  platform  since  both  of  these  genes  are 
directly  analyzed.  This  diagnostic  approach  is  very  effective  for 
broad-spectrum,  high-resolution  surveillance  situations.  The 
higher  per-assay  costs  of  the  RPM  are  offset  by  the  fact  that 
other  methods  would  involve  multiple  separate  tests  requiring 
more  time,  materials,  and  cost. 

The  use  of  two  additional  PCR-based  methods  illustrates  that 
other  assay  approaches  can  provide  similar  detection  of  avian 
influenza  virus  but  that  additional  tests  are  required  for  the  num¬ 
ber  of  pathogens  tested  for  by  RPM-Flu  to  be  matched.  While 
detections  of  the  presence  of  avian  influenza  virus  were  equiva¬ 
lent,  different  results  were  obtained  from  the  methods  for  sero- 
typing  11  of  the  37  samples.  Specific  PCR  remains  one  of  the  most 
inexpensive  molecular  diagnostic  methods,  making  it  useful  for 
identifying  the  presence  of  current  circulating  strains.  However, 
due  to  its  limited  scope,  it  is  not  easily  applied  for  identifying 
reassortments  or  new  variants  or  for  broad-spectrum  serotyping 
of  diverse  samples,  as  is  described  here  for  the  non-H5Nl  avian 
influenza  virus  samples  tested. 

RT-PCR-ESI-MS  requires  a  high  initial  capital  investment 
but  is  operationally  a  cost-effective  method  if  the  high-through- 
out  capabilities  of  the  platform  can  be  utilized.  Furthermore, 
this  approach  is  able  to  provide  more-detailed  information 
beyond  mere  identification,  such  as  its  capacity  to  track  the 
evolution  of  circulating  human  influenza  virus.  The  previously 
reported  capability  of  this  method  for  serotype  identification  is 
largely  confirmed  in  this  study  (30).  RT-PCR-ESI-MS  infer- 
entially  determines  the  identity  of  influenza  viruses  through 
base  composition  analysis  of  PCR  segments  of  polymerase 
(PB1,  PB2,  and  PA),  M,  NS1,  and  nucleoprotein  ( NP )  genes 
(30).  This  inferential  method  may  lead  to  the  discrepancies 
seen  in  samples  representing  less  prevalent  subtypes  with  lim¬ 
ited  sequencing  results  in  the  database.  Although  this  system  is 
capable  of  reporting  additional  information,  the  tests  and  anal¬ 
ysis  applied  to  the  samples  in  this  study  only  provide  serotyping 
information. 

While  the  RPM-Flu  chips  demonstrate  a  very  broad,  high- 
confidence  detection  and  analysis  capability,  the  technology 
can  be  improved  upon.  For  example,  the  current  arrays,  de¬ 
signed  to  provide  comprehensive  coverage  of  primarily  human 
bacterial  and  viral  pathogens,  may  not  effectively  provide  de¬ 
tailed  mutation  and  pathogenicity  information  for  all  subtypes 
of  influenza  virus.  This  impediment  reflects  the  limits  of  mi¬ 
croarray  real  estate,  which  demand  that  compromises  between 
the  depth  and  scope  of  coverage  must  be  made.  The  limitation 
will  become  less  of  a  factor  as  the  capacity  and  cost  of  rese¬ 
quencing  microarrays  improve.  In  addition,  no  conclusive  re¬ 
sults  could  be  obtained  from  the  critical  amino  acids  associated 
with  amantadine  drug  resistance  in  the  M2  protein.  This  is 


likely  due  to  the  fact  that  several  nucleotide  polymorphisms 
occur  within  the  same  18-base-pair  region  with  respect  to  the 
sequence  used  on  the  microarray.  More  specifically,  clustered 
single  nucleotide  polymorphisms  have  a  higher  probability  of 
not  being  identified  on  a  resequencing  microarray  (21).  There¬ 
fore,  only  ambiguous  calls  are  made  within  the  region,  and  no 
positive  evidence  of  the  mutation  is  provided. 

Future  iterations  of  RPM  design  can  alleviate  detection  of 
bases  in  clustered  mutation  hotspots  by  including  additional 
probes.  Another  issue  that  can  be  overcome  in  future  iterations 
by  modifying  design  strategies  and  analysis  algorithms  is  the 
current  difficulty  in  confirming  sequence  deletions.  Finally,  it 
should  be  noted  that  the  absolute  analytical  sensitivity  of  the 
human  influenza  virus  primers  is  slightly  higher  than  that  of  the 
avian  influenza  virus  primers,  suggesting  that  adjustments  in 
the  primers  may  improve  base  call  rates  by  providing  more 
amplicon.  To  maximize  the  benefits  of  using  a  resequencing 
microarray  for  avian  samples,  an  array  that  detects  many 
pathogens  of  concern  in  bird  samples  will  provide  more  infor¬ 
mation  and  could  be  placed  on  a  smaller,  less  costly  microar¬ 
ray.  Efforts  are  currently  ongoing  to  develop  such  an  array. 
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